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A new technique for the determination of the initial mass function in 
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ABSTRACT 

We present a new technique for the determination of the low-mass slope (ai; M* < 0.5M©) 
of the present day stellar mass function (PDMF) using the pixel space fitting of integrated 
light spectra. It can be used to constrain the initial mass function (IMF) of stellar systems with 
relaxation timescales exceeding the Hubble time and testing the IMF universality hypothesis. 
We provide two versions of the technique: (1) a fully unconstrained determination of the 
age, metallicity, and ot\ and (2) a constrained fitting by imposing the externally determined 
mass-to-light ratio of the stellar population. We have tested our approach by Monte-Carlo 
simulations using mock spectra and conclude that: (a) age, metallicity and ol\ can be precisely 
determined by applying the unconstrained version of the code to high signal-to-noise datasets 
(S/N=100, R = 7000 yield Aai w 0.1); (b) the M/L constraint significantly improves the 
precision and reduces the degeneracies, however its systematic errors will cause biased oi\ 
estimates; (c) standard Lick indices cannot constrain the PDMF because they miss most of 
the mass function sensitive spectral features; (d) the ot\ determination remains unaffected by 
the high-mass IMF shape (a^; M* J? 1M G ) variation for stellar systems older than 8 Gyr, 
while the intermediate-mass IMF slope (o^; 0.5 ^ M* < IM®) may introduce biases into 
the best-fitting ct\ values if it is different from the canonical value «2 = 2.3. We analysed 
observed intermediate resolution spectra of ultracompact dwarf galaxies with our technique 
and demonstrated its applicability to real data. 

Key words: methods: data analysis - globular clusters: general - galaxies: star clusters: 
general - galaxies: dwarf - galaxies: kinematics and dynamics - galaxies: stellar content 
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1 INTRODUCTION 

The stellar IMF is one of the most important intrinsic proper- 
ties of the star formation process. Its origin still remains unclear. 
Historically, it was thought to be a unimodal power law function 
dN/dM oc M~ {a ~^,a = 2.35 dSalpeterlfl955l) . The current 
working hypothesis (Kroupa 20Q2|) of the IMF universality sug- 
gests that its shape does not vary among different star forming re- 
gions and can be represented as a bimodal power law with differ- 
ent low- and high-mass slopes (ai = 1.3 for M < 0.5 Mq and 
o>2 — 2.3 at larger masses). Other possi bilities include variati ons 
of uni- or multi-moda l power laws (e .g. iMiller & Scalol [T979h or 
a lognormal function (Chabrier 2003). The IMF shape has a di- 
rect impact on galaxy evolution because it defines the number of 
Wolf-Rayet stars and supernovae and therefore the amount of star 
formation feedback. The low-mass IMF slope has strong influence 
on stellar mass-to-light ratios because of high M/L values of dwarf 

late type stars. 

Classical IMF determination techniques (see iKroupal <2002h 
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for a review) use direct star counts in open clusters and Hll asso- 
ciations. However, they are affected by various phenomena such 
as dynamical evolution of star clusters, variable and often strong 
dust extinction in star forming regions, as well as by the low num- 
ber statistics. In this light, compact stellar systems (CSS), that 
is globular clusters (G C) and ultracompact dwarf galaxies (UCD, 
Drink water et al.l 2003) which are orders of magnitude more mas- 
sive, free of the interstellar medium and populated by old stars 
present a unique laboratory to study the IMF. 

CSSs observed today might have experienced dynami- 
cal evolution effects on their stellar mass functions, i.e. the 
observed PDMF may differ from the I MF. It is known 
dSpitzerl [l987l; iBaumgardt & Makinj 120031; iKhalisi et alj l2007l : 
iKruiissen & Mieskej|2009f) that in globular clusters the dynamical 
evolution causes mass segregation, i.e. massive stars move towards 
the centre while low-mass stars migrate to the cluster outskirts, 
where they are tidally stripped during the passages close to the 
centre of a host galaxy or through its disc. This creates a deficit 
of low-mass stars in a cluster, changing the shape of its integrated 
stellar mass function. The characteristic timescale of this process 
is related to the dynamical relaxation time, which can be estimated 



2 Podorvanyuk, Chilingarian & Katkov 



M:ly 



Myr, 



for a CSS llMieske et alJ liijog) as t rcl = x - 

where A/d yn is in M© and the effective radius r e is in pc. Suffi- 
ciently massive GCs and UCDs have long relaxation timescales (a 
few Gyrs to many tnub) and therefore we can neglect the dynami- 
cal evolution effects in those systems. Hence, the PDMF there be- 
comes a good approximation of the IMF for low-mass stars which 
have not yet finished their evolution. 

One has to keep in mind that the mass segregation is not the 
only process which may result in the low-mass stars deficiency in 
CSSs. It can also be caused by the residual-gas expulsion from ini- 
tially mass segregated s tar clusters, as it was shown for galacti c 
GCs jMarks etalj|2008h and UCDs jDabringhausen et al.|[2oToh . 
However, the dynamical analysis of Fornax cluster UCDs does not 
show any sign s of the lack of low-mass stars in these relative ly mas- 
sive systems dMieske et alj|2008l : IChilingarian et alj|201 ll) . There 
is also a strong relationshi p between the low-mas s IMF slope and 
the concentration in GCs (De Marchi et al. 1 120071) . It turns out to 
be possible to estimate the phase of the low-mass star loss process 
from observations and therefore to make a conclusion about the 
equivalence of PDMF and IMF for every particular CSS. 

In this paper we describe a new technique of deriving the low- 
mass IMF slope by using the pixel fitting of intermediate- and high- 
resolution spectra of stellar systems integrated along the line of 
sight. This ap proach is inspired by the results of lChilin garian et al. 



2008, 20TTh where the stellar masses derived using the ISalpeterl 



19551) IMF turned to exceed the dynamical mass estimates result- 



ing in unphysical negative dark matter contents. 



2 NEW METHOD: SIMULATIONS AND APPLICATIONS 

We present two versions of the technique to derive the low-mass 
slope of the stellar IMF from integrated light spectra using full 
spectral fitting developed as an extension of the NBURSTS package 
jChilingarian et al.ll2007allbl) . 

NBURSTS is an approach to determine the parametrized line- 
of-sight velocity distribution (LOSVD) and star formation history 
(SFH) of unresolved stellar populations by means of the full spec- 
tral fitting in the pixel space. Previously the NBURSTS technique 
was mostly used with simple stellar population (SSP) models char- 
acterized by only two parameters, age and metallicity of the instan- 
taneous starburst event computed for a pre-defined IMF shape. In 
this study, we added an extra dimension to the grid of SSP models 
by varying the low-mass end slope ct\ of the Kroupa canonical IMF 
dKroupdl2002T) leaving the slope above 0.5A/q fixed (02 = 2.3). 
Hence, ai becomes an additional free parameter returned by the 
modified fitting procedure. We computed stellar population mod- 
els for "non-standard" ct\ values using both, high and low resolu- 
tion stellar atmosphere grids supplied with the PEGASE.HR evolu- 
tionary synthesis code (Le Borgne et al. 2 0041), the empiri cal high 
resolution ELODIE.3.1 stellar library jPrugniel et ail2 007 ) and the 
low resolution BaSeL theoretical atmos pheres (ILejeune et alj 19971) 
used in the original PEGASE.2 package I IFioc & Rocca-Volmerangej 
1 1997b . The full spectral fitting is performed with high-resolution 
models while the mass-to-light ratios of stellar populations in broad 
band filters are estimated using the low-resolution grid. We used 
standard PEGASE.HR settings for all other parameters: 5 per cent of 
close binary systems, low and high cut-off stellar masses of 0. 1 and 
120 Mq correspondingly. 

In the first version of our new approach an observed spectrum 
is fitted against an optimal template represented by a linear combi- 



nation of SSPs each of them characterised by the age, metallicity 
and the low-mass IMF slope ai, determined in the same minimiza- 
tion loop, using the 3-dimensional cubic spline interpolation on the 
pre-computed grid of high-resolution SSP models. 

In this paper we deal with observed spectra of CSSs which 
can be well represented by single component SSP models, how- 
ever the code provides a possibility of fitting multiple component 
stellar populations. The model grid has 25 nodes in age (10 Myr 
to 20 Gyr), 10 nodes in metallicity ([Fe/H] from -2.5 to 1.0 dex) 
and 10 nodes in the low-mass IMF slope (0.4 to 3.1 with a step of 
0.3). During the fitting procedure, th e models are broadened using 
the Gauss-Hermite parametrization dvan der Marel & Franx|[l993f) 
of the LOSVD which is penalized toward purely Gaussian solu- 
tion as proposed by ICappellari & Emselleml J2004h . In this paper 
we consider only purely Gaussian LOSVDs. 

The x 2 value (without penalization) is computed as follows: 



2 (F l -Pi P (T i {t,Z,a 1 )®jC(v,o-,h 3 ,h i ) + P 2q )) 
X = > , (!) 



AFi- 



where Fi and AFi are observed flux and its uncertainty; C is 
LOSVD used as a convolution kernel for Ti(t, Z,ai), the flux from 
a synthetic spectrum, represented by a linear combination of SSP 
models characterised by their ages t, metallicities Z, and low-mass 
IMF slopes oti, and convolved according to the line-spread func- 
tion of the spectrograph; Pi p and Pi q are multiplicative and ad- 
ditive Legendre polynomials of orders p and q for correcting the 
continuum; v, a, /13 and /14 are radial velocity, velocity dispersion 
and Gauss-Hermite coefficients of the LOSVD respectively. 

The x 2 minimization is performed using the constrained non- 
linear Levenberg-Marquardt minimization implemented in the MP- 
FIT IDL package (by C. Markwardt, NASA). Since it requires the 
second derivatives to be continuous, we use a three-dimensional 
spline interpolation of the SSP grid to evaluate a model for a given 
set of (t, [Fe/H], «i). In this fashion, the values of t, [Fe/H], and qi 
are returned by the minimization procedure in the same loop along 
with the parameters of the continua. 

The second version of our technique is mathematically iden- 
tical to the pure NBURSTS technique but the specific grid of input 
SSP models is constructed and supplied for every observed spec- 
trum. Here we assume that the object being studied contains no 
dark matter. This assumption is probably true for GCs and UCDs. 
In this case, the dynamical mass-to-light ratio (M/L)d yn equals to 
the stellar mass-to-light ratio (M/L),. For a number of compact 
stellar systems, the (M/L)a yn values were determined from the 
analysis of their internal structure and observed velocity dispersion 
profiles and are available in the literature. On the other hand, we 
derive the (M/L), values from low-resolution stellar population 
models for every set of (t, [Fe/H], Qi) and for every given t and 
[Fe/H] in the old stellar population regime this function is mono- 
tonic. Hence, if we know (A// £)d yn and impose the zero dark mat- 
ter constraint, for every point on the (t, [Fe/H]) plane we can find 
the value of ai such as (M/L) * in that point of the parameter space 
equals to (M/L)d yn - Thus, we can compute a grid of SSP models 
in the age-metallicity space varying the low-mass IMF slope so that 
the (Af/L)* values are constant all over the grid and equal to the 
"external" (A//L)d yn . Along with this grid, we will also map the 
behaviour of et\ as a function of t and [Fe/H]. Then, if we feed this 
SSP grid to the standard NBURSTS full spectral fitting procedure 
and determine the pair of best-fitting values of age and metallicity, 
we will automatically get the ai value corresponding to this best- 
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fitting solution which will measure the low-mass PDMF slope of 
the stellar population being studied. 

It may happen that in order to reach the imposed value of 
(M/L) the low-mass IMF slope has to be outside our model grid 
(0.4 < ai < 3.1). In this case, we take the model corresponding 
to the nearest limiting ot\ value (maximal or minimal). 

In some cases, e.g. high-resolution spectra of extragalac- 
tic GCs or UCDs where the internal velocity dispersion can be 
measured from a spectrum itself and where it uniquely defines 
(A//L)d yn , it is possible to modify the minimization procedure 
and close the loop so that no external (M/L) value needs to be 
supplied. However, in this case the procedure becomes very com- 
plex and we do not consider this approach in this paper. As we 
demonstrate from simulations, velocity dispersion does not exhibit 
any degeneracies with age and oti, therefore if a is determined in- 
dependently, it should not affect the final ai estimate. 



2.1 Full spectral fitting of mock data 

We performed Monte-Carlo simulations in order to study the pre- 
cision, stability, and possible biases of the proposed technique for 
the determination of the low-mass IMF slope. We generated mock 
datasets (10 to 1000 noise realisations for every combination of in- 
put parameters) for the signal-to-noise ratios of 25, 50, and 100, 
and the spectral resolution values R — 4000 . . . 10000. Then we 
fitted them using both versions of our technique. 

In Fig. Q] we present the results of the full spectral fitting 
of 1800 mock spectra (300 noise realisations for the ages t=8, 
10 and 12 Gyr; Z = —0.3 dex; ai =1.3 and 2.2, i.e. close to 
Kroupa and Salpeter IMFs) for the spectral resolution R = 7000 
and the signal-to-noise ratio S/N=100 in the wavelength range 
3900. . .6800 A. The results are presented in a way showing the de- 
generacies between derived stellar population parameters and the 
qi IMF slopes. The results of the simulations made using metal- 
poor stellar populations ([Fe/H]=— 1.0 dex) look very similar ex- 
cept the ~30 per cent higher uncertainties. 

Here we see that even the first version of the technique is 
able to recover the Qi values with the uncertainty of about 0.1- 
0.2. When we impose the M/L ratio, the precision significantly 
improves and the uncertainties of ai become better than 0.1. We 
see the strong age-ai degeneracy because the higher mass fraction 
of low-mass stars in older populations may compensate the lower 
derived values of ai. The accuracy of about 0.06 well compares 
to the best available IMF determination for Galactic open clusters 
Pleiad es (Aqi=0.15) and M35 (Aqi=0.12) presented in lKroupal 
d2002h . 

In order to test the sensitivity of the second version of our tech- 
nique to the biased input M/L values, we performed the simula- 
tions increasing and decreasing the M/L expected from the stellar 
population models by 10 per cent. The results are shown in Fig.Q] 
by green (0.9 ■ M/L) and cyan (1.1 • M/L) contours. It demon- 
strates how M/L systematic errors translate into ai biases. One 
has to keep in mind that usually the dynamical mass of a real stel- 
lar systems is known with the precision worse than 10 per cent. 

2.2 The analysis of the distribution of the IMF sensitive 
information 

For the practical usage of our technique and future impovements 
of the observational strategy to study the IMF in real stellar sys- 
tems, one needs to know which spectral features have the highest 



sensitivity to the IMF shape. In order to identify them, we analysed 
the distribution of the ai -sensitive information in 20A-wide bins 
in the optical spectral domain for metal-rich and met al-poor stel- 
lar pop ulation in a similar fashion to that proposed in lChilin garian 
(2009) for the analysis of age and metallicity sensitivity. For our 
test (see Fig. [2] we chose two high resolution (7? = 10000) PE- 
GASE.HR SSP models in the wavelength range 3920. . .6700 A with 
the age t — 10 Gyr and ai=1.3 having metallicities [Fe/H]=— 1.0 
and —0.3 dex. Then, we convolved them with Gaussian kernels 
corresponding to the internal velocity dispersions of 10 km s _1 
and thus obtained two model spectra or "reference SSPs". Then for 
every of them we varied the ai value by 0.1. Later, we fitted these 
models against their "referenc e SSPs" using the pPXF procedure 
(Cappellari & Emselleml l2004l) with the 10th order multiplicative 
polynomial continuum. The total amount of "information" equals 
to the sum of partial derivatives of \ 2 as a function of Qi over all 
pixels normalised to 100 per cent. The Nal D line region shown in 
yellow in Fig.[2]is excluded from consideration because it is a fea- 
ture in stellar spectra heavily affected by the interstellar absorption 
line and therefore cannot be reliably modelled in the evolutionary 
synthesis. 

We see that the distributions are very similar for metal-rich 
and metal-poor populations (red and blue histograms respectively). 
An exception is the Ha line (6563 A) which contains considerable 
amount of the information for metal-poor populations and is less 
significant at higher metallicities. Most of the IMF-sensitive spec- 
tral regions are those containing numerous but relatively faint ab- 
sorption lines sensitive to the surface gravity in the atmospheres of 
late type (GKM) stars. These lines are concentrated on both sides 
of H/3 (4760-4800A and 4900-4960A) and between 5400A and 
6400A with the most sensitive part around 6100A associated with 
Cal, Vl, Col, Nil, and Til absorption lines. Worth mentioning that 
when varying log g, the Cal lines change in the opposite direction 
compared to all other elements from the list above. Details will be 
given in a separate paper. 

2.3 Effects of the high-mass IMF shape on the ai 
determination 

In this study we use the Kroupa canonical IMF tooupall200l) bv 
varying the low-mass end sl ope ct\ and le aving the high-mass slope 
fixed (q 2 = 2.3). Recently, Mar ks et al.l d2012h have demonstrated 
that stellar systems which underwent high star formation volume 
densities at a parsec scale require the introduction of 03 for massive 
stars (> 1A/q). Using Monte-Carlo simulations we have checked 
how changes of otz influence the fitting results for a "standard" bi- 
modal canonical IMF. We generated 4500 spectra (300 noise reali- 
sations for the ages t=4, 6, 8, 10 and 12 Gyr; as=l. 7, 2.0 and 2.6, 
ai=1.3, [Fe/H] — —1.0 dex). Then we fitted them using the first 
version of our technique with a model grid having 03 = 2.3. The 
results of this fitting are presented in Fig.|3]as data points connected 
by solid lines. 

Here we clearly see that the Qi determination remains un- 
affected by 03 values for stellar systems older than 8 Gyr. It is 
trivially explained by the fact that the lifetime of massive stars 
(> IMq) is shorter than 8 Gyr therefore they do not contribute 
to the integrated light of a stellar system. The conclusion from this 
experiment is that any stellar system of the age of 8 Gyr or older 
can be analysed with our standard set of the SSP models with the 
fixed value of 03 = 2.3 when using the first version of our IMF 
determination technique. 

We performed similar analysis in order to test how the best- 
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Figure 1. Results of the full spectral fitting of mock spectra. Red and blue contours correspond to the first (all free parameters) and the second (constrained 
M/L) versions of the technique. Cyan and green contours demonstrate how biased values of M/L ratios affect the results in the constrained version of the 
technique. Contours correspond to the 1-, 2- and 3-cr levels. 
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Figure 2. The distribution of a\ -sensitive information along the wavelength in 20 A-wide bins in the high-resolution (R=70()0) spectra of metal-rich (red, 
[Fe/H]=-0.3 dex) and metal-poor (blue, [Fe/H]=-1.0 dex) stellar populations. The distributions are normalised to the unity (100 per cent). The stellar population 
spectra are shown above. The Nal D region shown in yellow is excluded from consideration (see the text). Green areas denote the Lick index definitions. We 
clearly see that apart from a few exceptions, the IMF-sensitive information is not associated with Lick indices (see the text). 



IMF determination in unresolved stellar populations 5 



2.5 



1.5 



a =1.7 




Figure 3. The ct\ determination biases introduced by the variations of a.2 
and 03. For every set of parameters, error bars show the standard deviation 
of the best-fitting a± calculated from 300 noise realisations in every point. 
Solid and dash dot lines are for the fitting results when varying 03 and a.2 
respectively. Arrows indicate that the lower limit of the model grid (ai = 
0.4) was reached for some mock spectra. 



fitting a\ values are affected by the variations of a%. We generated 
3000 mock spectra (300 noise realisations for the ages t=4, 6, 8, 10 
and 12 Gyr; a2=2.0 and 2.6, Qi=1.3, metallicity Z — —1.0 dex). 
Then we fitted them in a similar fashion to the Q3 test. These results 
are also shown in Fig.|3]as data points connected by dash dot lines. 

We see that 02 variations introduce significant biases into 
the best-fitting ai values. The bias is larger for younger popula- 
tions (ai=2.0 instead of 1.3 for 4 Gyr) than for older populations 
(ai=1.7 instead of 1.3 for 12 Gyr). This fact is also easy to explain 
by the theory of stellar evolution: for older stellar populations all 
massive stars have already evolved into the remnants and the con- 
tribution of stars with masses M > 0.5Mq decreases as the age 
increases. The lower a.2 values correspond to the excess of these 
stars compared to the canonical IMF shape and when we fit such 
a spectrum with the models computed for a.2 = 2.3, it requires 
shallower low-mass end slope a\ to partially compensate the lack 
of more massive stars which were present in the model because of 
lower «2 • 

However, we stress that according to recent studies there is 
no evidence for Q2 to vary: "The IMF of massive stars is well de- 
scribed by a Salpeter/Massey slope, 02 = 2.3, independent of envi- 
ronment as deduced from resolved stellar populations in the Local 
Group of galaxies. Unresolved multiple stars do not si gnificantly 
affect the massive-star power-law index" (Krou paet al.l201 lh . 



2.4 Full spectral fitting of real data 

We applied our technique to observed spectra having different 
signal-to-noise ratios, spectral resolutions and wavelength cover- 
ages. We used three publicly available datasets for compact stellar 
systems : (1) intermediate re solution (R=1300) Galactic GC spec- 
tra from lSchiavon et al.l a2005l> obtained with the R-C spectrograph 
at the 4 m Blanco telescope at CTIO in the wavelength range 
3900< A <6500 A; (2) high-resolution (R=17000) intermediate 
signal-to-noise spectra of UCDs in the Fornax cluster obtained with 
FLAMES/Giraffe at ESO VLT in the HR9 se tup covering a nar- 
row wavelength interva l, 5120< A <5450A jMieske et al]|2008l ; 
IChilingarian et al.l201 lh downgraded to R=10000, and (3) interme- 
diate resolution (R= 1500-2000) high signal-to-noise ratio spectra 
(3900< A <6800 A) of Fornax a nd Virgo cluster UC Ds obtained 
with Gemini GMOS-N/GMOS-S ^Francis et alj2012h . 



The dynamical M/L measurements fo r some of 
ISchiavon et al.l J2005t) spectr a were available from iDubath et al.l 
dl997h and lLane et all fcoiOh . However, none of the clusters had 
^rel ^ ^ Gyr, and those spectra were integrated and extracted only 
in the cluster core region. In addition, classical long-slit spectra 
of GCs used for the dynamical analysis and consequently our MF 
analysis suffer from systematic uncertainties caused by the limited 
number of individu al stars falling into the slit (see discussion in 
Dubat het alJI 19971) . Another problem is that the nuclear regions 
will have higher weights over peripheral parts because the cluster 
light fraction falling into the slit decreases as R~ . Therefore, we 
could not obtain any reliable estimates of qi values in these objects 
but only use them to evaluate statistical errors of the techniques 
for spectra of such a type. The statistical uncertainties of the first 
version of our technique turned to be about Aqi ~0.3-0.5. Fixing 
the M/L values reduced them down to 0. 15. 

Spectra of extragalactic clusters and UCDs have a number of 
advantages over Galactic GC spectra although they are more dif- 
ficult to obtain because the targets are much fainter. The main ad- 
vantage is that these spectra will be very close to the integrated 
spectra because at the distance of Virgo and Fornax CSSs usually 
stay unresolved from the ground. Therefore, there will be no sys- 
tematic errors connected to the low number statistics of stars on the 
slit, and also there will be no R^ 1 weighting of peripheral parts. 
Another advantage is that in rich environments such as Fornax and 
Virgo, the number of GCs is so high that they populate very well 
even the brighter end of the luminosity function, therefore it be- 
comes easy to find objects with t rc i ^ tuub which did not undergo 
the dynamical evolution and therefore are suitable for the determi- 
nation of IMF and not only the PDMF 

Our experiments with the high-resolution FLAMES/Giraffe 
spectra of Fornax cluster GC/UCDs demonstrate that those data can 
be used for the PDMF determination only if the dynamical M/L 
values are known from the analysis of their internal dynamics and 
structural properties. Then the statistical uncertainties of the M/L 
constrained fitting stay in the range of 0.05. . .0.1 depending on the 
S/N ratio and the overall errors are systematic and fully imposed 
by the M/L uncertainties. However, very short wavelength range 
not including any strong IMF sensitive features makes the ai de- 
termination impossible for the first unconstrained version of our 
technique bringing the statistical error to Aqi =0.8 ev en for the 
brightest object, UCD 3 (F 19 in lChilingarian et al.ll201ll) . 

Analysis of high S/N integrated spectra for Virgo/Fornax 
UCDs obtained with Gemini GMOS-N and GMOS-S spectro- 
graphs in a broad wavelength range demonstrates a reasonably 
good agreement between the two versions of our technique (see 
Fig EJ- We used published dynamical M/L ratios for two V irgo 
cluster UCDs (VUCD 3 and VUCD 5, Evs tigneeva et al.l2007l) and 
two F ornax cluster UCDs (UCD 1 and UCD sjchilingaria n et al.l 
|2008J). The adopted (M/L) y( jyn values in the Solar units are: 
5.4 ± 0.9, 4.1 ± 0.8, 5.0 ± 0.6, and 3.4 ± 0.9 for VUCD 3, 
VUCD 5, UCD 1, and UCD 5 correspondingly. For UCD 1 and 
VUCD 5 both, fully unconstrained fitting and the Af/L-imposed 
version of the technique agree well within the statistical uncertain- 
ties of the unconstrained fitting. For UCD 5 having low metallicity 
and therefore much weaker absorption lines, the unconstrained fit- 
ting results in very large uncertainties while the adopted (M/L)v 
value puts it close to the edge of the model grid on the MF slope 
(2.8). For VUCD 3 the situation is the opposite, the unconstrained 
fitting yields ai ~ 2.9, however one has to keep in mind that this 
object has extremely h igh value of the [a/Fe] chemical abundance 
ratio dFrancis et al.ll2012f) resulting in the significant template mis- 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
[Fe/H], dex a 

Figure 4. Determination of the low-mass IMF slope using both versions of our technique for four UCDs: UCD 1 and UCD 5 from the Fornax cluster; VUCD 3 
and VUCD 5 from the Virgo cluster. Blue and red crosses and ellipses denote the best fitting values and their 1-cr statistical uncertainties computed from 
the covariance matrix for the unconstrained and M/ L constrained versions of our technique correspondingly. Red boxes demonstrate the range of systematic 
errors for the constrained fitting corresponding to i-cr uncertainties of the adopted M/ L ratios. Black lines on the lower panel show the limits of our grid of 
SSP models. 



match, therefore the unconstrained measurements may get biased. 
At present, we do not have models required to study the biases 
which may be introduced to the MF measurements from integrated 
light spectra by significantly non-solar abundance ratios. 



3 DISCUSSION AND SUMMARY 

We presented a novel technique for the determination of the low- 
mass end slope of the PDMF of unresolved stellar populations us- 
ing pixel fitting of spectra integrated along the line of sight. In stel- 
lar systems with long dynamical relaxation timescales such as mas- 
sive GCs, UCDs and galaxies, where one can neglect the effects of 
the dynamical evolution, the PDMF is identical to the IMF at low 
masses. Our technique is available in two versions having different 
applications: (a) fully unconstrained determination of age, metal- 
licity and the low-end PDMF slope ai which can be used to anal- 
yse high signal-to-noise spectra of CSSs and normal galaxies given 
that the wavelength coverage includes IMF sensitive features; (b) 
constrained determination of the same parameters by imposing the 
external value of the M / L ratio which can be, for example ob- 
tained from the dynamical modelling and therefore assuming the 
zero dark matter content in a stellar system being studied. 

We used Monte-Carlo simulations with mock spectra to test 
both approaches and conclude that: (1) age, metallicity, and IMF 
can be precisely determined in the first unconstrained version of 
the code for high signal-to-noise datasets (e.g. S/N=100, R=7000 
give the uncertainty of at\ of about 0.1-0.2); (2) adding the M/L 
information significantly improves the precision and reduces the 
degeneracies, however systematic errors in the imposed M/L ratio 
will translate into offsets in ai . 



We performed the analysis of the IMF-sensitive information 
distribution across the optical spectral domain. One can clearly see 
in Fig. [2] that apart from a few exceptions, IMF sensitive informa- 
tion is not associated with prominent absorption line features (and 
line strength indices corresponding to them, e.g. Lick indices, see 
Worth evf 19941) where a significant fraction of age- and metallicity- 
sensitive information is contained. The exceptions are: the Cal fea- 
ture (A = 4227 A) and the Mg6 triplet (A « 5172 A). However, 
these two indices are very sensitive to the metallicity and the [a/Fe] 
abundance ratios and therefore can hardly be used to constrain the 
IMF. Another IMF sensitive Lick index is Fe5782 which is known 
to be difficult to measure and calibrate. The only strong IMF feature 
left is the Ti02 index including about 7 per cent of the total IMF- 
sensitive information. Hence, "standard" absorption line strength 
indicators cannot be used to determine an. 

According to the recent studies dMarks et al.ll2012h , the stel- 
lar IMF at high masses (> IMq) appears to depend on the 
star-formation rate density on a parsec scale and even becomes 
top-heavy when it surpasses O.IM© / (yr • pc 3 ). It may be- 
come increasingly bottom-heavy at high metallicity. The increased 
number of compact stellar remnants such as neutron stars and 
black holes for a top-heav y IMF will elevate the M/L values 
jPabringhausen et al.l I2009D and mimic the dark matter in CSSs. 
Conversely, for a top-light IMF the M/L will be lower then for 
the canonical Kroupa IMF. This will become a caveat for the sec- 
ond (M/L constrained) version of our technique because the M/L 
ratios computed from the stellar population models based on the 
Kroupa canonical IMF will be incorrect. We quantified this effect 
by comparing the PEGASE.2 predictions of stellar M/L ratios for 
stellar populations having different tvi, as, and ages. As expected, 
the a3-induced M/L variations are the most important for bottom- 
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light populations (qi < 1) because of the higher relative mass frac- 
tion of stellar remnants. However, the overall effect of the 03 vari- 
ations on M/ L ratios is lower than that of qi for the same amount 
of the IMF slope change. For instance, going from as = 2.3 to 
Q3 = 1.7 changes the V band M/L ratio by 11 (for ct\ = 3.1) to 
26 per cent (for a\ = 0.4) for ages exceeding 10 Gyr. For 8 Gyr 
the effect is even 2 per cent smaller. The top-light IMF (03 = 2.6) 
lowers the M/L by as little as 1 to 6 per cent for ages above 8 Gyr. 

At the same time, the unconstrained version of our fitting tech- 
nique will provide correct estimates of the low-mass IMF slope 
if applied to older stellar populations (t > 8 Gyr) where mas- 
sive stars affecting the variable high-end IMF slope have already 
evolved into stellar remnants (see section 2.3). It also can be ap- 
plied to the data for galaxies because the assumed dark matter con- 
tent does not have to be zero. The unconstrained version can be 
applied to check whether the IMF is bottom heavy in massive ellip- 
tical galaxies wh ere low-mass stars might have efficiently formed 
in cooling flows dKrou pa & Gil morelll994l) . However, we have to 
keep in mind that the difference between SSPs and real star forma- 
tion histories of galaxies as well as significantly non-solar element 
abundance ratios may introduce some biases. 

We applied our technique to observed spectra of CSSs and 
demonstrated that even for low resolution datasets (R = 1500) 
with high S/N-ratios the IMF slope can be precisely determined 
even using the unconstrained fitting of ai , t , [Fe/H] . The M/L con- 
strained approach provides even higher precision and agrees well 
with the unconstrained version for objects without peculiarities in 
stellar populations and chemical abundances. 

High spectral resolution will improve the precision only if 
IMF-sensitive spectral features are covered. Extragalactic GC sys- 
tems where structural parameters of clusters are available from the 
Hubble Space Telescope data will have advantage over the Galac- 
tic GCs because the integrated velocity dispersion measurements of 
such systems are much easier to obtain for large samples of targets 
using multi-object spectroscopy. 

We conclude that by applying our technique to high-quality 
optical observations of CSSs we are able to reach better precision 
of the IMF determination than that made with direct star counts in 
nearby open clusters and check the IMF universality hypothesis. 
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